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We analyze the conditions that lead to an adiabatic evolution with a geometric phase accumulated. To real- 
ize the adiabatic evolution, the instantaneous eigenenergies and the degeneracy of the Hamiltonian can change 
rapidly during the evolution. Using pulse sequences or fast modulation fields, we can suppress the state transi- 
tions and realize the adiabatic evolution in a much shorter time. 
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Introduction. The adiabatic theorem in quantum mechan- 
ics concerns the evolution of quantum systems subject to 
slowly varying Hamiltonians |1|. It says that the transitions 
between the instantaneous eigenstates of a Hamiltonian is 
negligible if the change of the Hamiltonian is much slower 
than the energy gaps between the instantaneous eigenstates. 
Berry found that in the adiabatic evolution there is a geo- 
metric phase determined only by the path in addition to the 
dynamic phase [2]. Wilczek and Zee generalized the result 
to non-Abelian geometric phase for degenerate Hamiltoni- 
ans 01 . The geometric phase is useful in quantum informa- 
tion processing |4-6|. When the change of the Hamiltonian 
is not much smaller than the energy gaps, the transitions be- 
tween the eigenstates conceal the geometric phase, that is, 
these transitions induce an error evolution U^r on the quan- 
tum system. Obtaining an exact expression of the error evolu- 
tion £/En- is important. By Ueit, we can quantify the accuracy 
of the adiabatic approximation and solve some recent debates 
(e.g., Refs. II7 UT21 ') on the conditions and validity of adiabatic 
theorem. 

Dynamical decoupling is an efficient method in suppress- 
ing the error evolutions including decoherence and unwanted 
state transitions, through fast operations on the quantum sys- 
tem I1T3T - IT71 . The performance of dynamical decoupling was 
demonstrated by experiments (e.g., Refs. [18 -22]). It would 
be of practical interest if we can use the idea of dynamical 
decoupling to realize adiabatic evolution in a shorter time by 
suppressing undesired transitions. 

In this Letter, we develop a general formalism to quantify 
the error evolution Ueu of the quantum system in the adia- 
batic evolution. In the adiabatic limit, the error evolution ap- 
proaches an identity operator (1/rjrr — > !)• From the expres- 
sion of (7eit, we find that to realize Ueu — > L the Hamiltonian 
can change rapidly and the energy levels of the instantaneous 
eigenstates can have degeneracy or energy crossings during 
the evolution. Using dynamical decoupling control, we can 
suppress the state transitions and realize the adiabatic evolu- 
tion in a much shorter time. Given a path of a parameter space 
for adiabatic evolution, we can realize the adiabatic evolution 
by applying a sequence of control pulses which sample the 
path for adiabatic evolution (see Fig. [T] for an illustration). 
The time intervals between the pulses can be small, and the 
adiabatic evolution can be achieved in a very short time. 




FIG. 1: (Color online). A sequence of pulses applied to a spin- ^. The 
red dots indicate the directions of the pulses on the Bloch sphere. The 
blue semi-circles illustrate the directions of spin rotation. The green 
dashed line shows the effective path of an adiabatic evolution. This 
successive rotations induce a geometric phase. 

General formalism. For simplicity we consider quantum 
systems of finite dimension, and our formalism can be gener- 
alized to the systems with continuous spectra. The quantum 
system is subjected to a time-dependent Hamiltonian H(t) = 
H(R), where the parameters R = R(t) = (R l (t),R 2 (t), ■ ■ ■ ) 
change in time as 

dR = <xidt. (1) 

The evolution of arbitrary quantum states (including mixed 
states) is described by the evolution operator U(t) = U{R), 
which satisfies the Schrodinger equation (with h — 1) 

Au(f) = iU(t) = H{i)U{t), (2) 
dt 

where a dot denotes a time derivative. We will show 
that at the moment T the evolution operator U(T) = 
UDyn(T)UGeo(T)UEn:{T) is the product of a dynamic phase 
factor Uu yn (T), a geometric phase factor UG eo (T) = UceoiRr) 
with Rt = R(T), and an error evolution Ueit(T). For an adia- 
batic evolution, Ubb(T) — * I- 

The instantaneous orthonormal eigenstates \n'.) = \n At)) = 
|«*) of the Hamiltonian at time t satisfy 

H(f)\n'j) = E nJ (f)\n )) = E nJ (R)\n* >. (3) 
The subscript j may be used to label degenerate states. 
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Substituting the transformation U(t) = S(t)Us(t) in Eq. (|2| 
with the unitary operator S (t), we obtain 

ills (f) = S '"(f) [H(t) - iS (t)S *(*)] 5 (t) tf s (f). (4) 

Let 5(f) = Ux) yn {t)Uc\(.t), and we extract the dynamic phase 
factor of the system evolution by 



cM')=Z e "'- c ' ! " j " > * | "j >< "; 1 
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j JK o \nj)(nj\ 



(5) 



And the unitary operator 

f/Gl(0 = f/Gl(/?) = 2K>< n J°l 



(6) 



depends only on the geometric property of the eigenstates. We 
have used the notation Ro = R(0). The solution of Eq. Q at 

time T is U S (T) = Te^ti ' Hs( - t)dt with H s = H G2 + H Err , where 
T is the time-ordering operator and 

HmV) = - J] InlKn^iny^-^, (7a) 
",p,q 

H E u(t) = - KHn^Kmy^-^. (7b) 
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We write {/$ in the interaction picture of H G2 as Us(t) - 
U G2 (t)U En (t), where t/ G2 (f) = Te^***™* and t/ ElT (f) = 
Te -i£H En (s)ds wkh 5firW = ul 2 (t)H ETr (t)U G2 (t). When the 

effect of the oscillating factors e l h { - E "- p ~ E "-' l)ds in Eq. |7a| van- 
ishes, H G2 (t) becomes geometric. One way to satisfy this con- 
dition is to make the energies E„ )P (t) = E nq (f) the same for all 
state labels n, p and q, which is just the usual consideration 
of degenerate Hamiltonians Q. We note that under a more 
general condition 



Jo 



(£„, p - E n , q )ds = 0, +2n, ±4n, 



(8) 



during the system evolution, we still have e l ^ Ea - p ~ En -^ ds = 
1. And this is possible when we use pulse sequences 
and the dynamic phases J E„ p (t)dt change suddenly. With 
e > f (E„.i,-E„. q )ds _ ^ ^ p ro p a g a tor U G2 (t) is a geometric phase 
factor determined by the parameter R at time t, 



U C2 (R) = Pe 



i Jg* £*, M \A ° ><»f l«f )(>"" W 



(9) 



where V is the path-ordering operator and Vr< = 
(^|r, • ■ ■ ) acting on |n* '). Then the whole evolution op- 
erator 



= U Dyn (T)U Geo (T)U En (T), 



is the product of the dynamic phase factor Ur> ya (T) in Eq. d5j, 
the geometric phase factor 



f/ Ge o(7') = U Geo (Rr) = U gi (Rt)U G2 (Rt) 



and the error evolution 



in 



(11) 



U BtT (T)=Pexp 



r K T 
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(R)H n %(K) ■ dR 



(12) 



where 



H p n %(R) = ^ 2 (if)|n°><n*]iV«|mf><m°]f/ G2 CJf) 

= ut(R)\n*){n*\H R \m«)(m«\U Geo (R), 



(13) 



depends on the geometry of the path and the modulation func- 
tion 

The operators H P %(R) describe the effects of the transitions 
\ m q ) l n p) between different subspaces. When the fre- 
quencies of the Fourier components of the modulation func- 
tions F p 'f n (R) are much larger than the typical frequencies of 
Hf t %(R), the effects of the transitions \m^) <-» |«*) are sup- 
pressed and U En {T) — > /. 

One advantage of writing the whole system propagator 



U(T) in the form of Eq. ( 10 1 is that the error in the adiabatic 
approximation is completely described by the error evolution 
U En -(T). 

In this Letter, we measure the error evolution by the average 
deviation 



au Eii - = \mu En (T) - 



(15) 



where the over bar is the average over all possible initial states 
AU En - -» when U En (T) -» /. 
If we are interested in the adiabatic evolution in a sub- 
space, we introduce the projection operator II on the asso- 
ciated eigenspace and the projection operator out of the sub- 
space E = / - EL with n 2 = II and E 2 = E. Here the asso- 
ciated eigenspace of E may contain continuous spectra when 
our formalism is generalized to the systems of infinite dimen- 
sion. The adiabatic evolution in this subspace is realized when 

Uarirya -> n. 

Adiabatic evolution. Let the eigenvalues of the Hamilto- 
nian E, hp (t) - E n (t) for all n and p, one case of the condition 
Eq. ([SJ. The usual adiabatic condition requires that the change 
of the set of parameters R is very slow compared with the en- 
ergy gaps E„ - E m t 0, i.e., - E m )\ = \oj/(E„ - E m )\ 
very small. The slowly varying Hamiltonian implies that the 
dynamic phase factors F P %(R) [Eq. ( 14 1] in Eq. (TT2b are fast- 



oscillating functions. Note that the operators H n ' ln (R) are ge- 
ometric and they vary slowly during the evolution. There- 
fore the fast-oscillating phase factors F p 'Jj,(R) make the con- 
(10) tribution from H p f n (R) negligible. In this adiabatic limit 
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\(of(E„ - E m )\ — > (i.e., (x> — > with finite energy gaps 
\E n - E m \ + 0, or finite a> + with |.E„ - £„,| — > oo), we 
have Ueu:(T) — > I. This is the statement of the conventional 
adiabatic theorem. 

Geometric phase revealed by dynamical decoupling. The 
expression of U(T) by Eq. ( 10 1 gives us the following obser- 
vation: if we somehow realize U^ n (T) — > / with the geomet- 
ric phase UGeo{T) the same as that in an adiabatic evolution, 
we will realize that adiabatic evolution. Note that the key to 
make U^iT) — » / is the existence of the modulation function 
F^f n (R) in Eq. ( 12 1. The general condition to realize adiabatic 
evolution is that F l ^f n (R) is a modulation function that makes 



the effect of H%f„(R) in Eq. ( 12 1 vanish for n + m. We will 



show later that, energy degeneracy of E np and E„ uq is actually 
allowed during the system evolution. This is quite different 
from the conventional belief that we need E n>p + E m<q during 
the adiabatic evolution when m + n. 

In the theory of dynamical decoupling, there are similar 
modulation functions to make the error evolution vanish 11231 . 
This observation leads to the realization of adiabatic evolu- 
tion by dynamical decoupling control. To illustrate the idea, 
we consider a quantum system controlled by a sequence of N 
unitary pulses, 



m ^\n R "){n R " 
1 j ' v ; 



(16) 



»■./ 



during some time T, which can be arbitrary short. The pulses 
are applied in the order of the parameters R^ = Ri, R2, ■ ■■ , 
Rjy. For simplicity, we first assume that between the pulses 
the system Hamiltonian H(t) = 0. To have UGi(t) being com- 
pletely geometric, we require the condition Eq. (jH}, i.e., 

Yj KiW - 0n, q (Rv)] = 0, +2tt, +47r, . . . , (17) 

for all n, p, q, and A. One solution is 6 np (Ru) = 6 n q {R^). 
If we want to have M non-degenerate subspaces, we may pe- 
riodically apply pulses with GnjiR^ - 2nn/M 112411251 . An 
example of the pulse is the rotation of an angle 2n/(2J +1) 
on a spin-/ by a magnetic field with a constant direction that 
defines the eigenstates !«*'')■ 

The dynamic phase is ^ OnjiR^) and the geometric phase 
factor f/ceo(^r) is given by Eq. (Hi with the path sampled by 
the points R^. In the error evolution UBafX), the modulation 
function becomes F^%{R) = g' ^K/W-fW*,,)^ wnere ^ is 
the summation over the phases induced by the pulses that have 
been applied. The aim of the pulse sequence is to suppress the 
effect of H%%(R). 

Note that for a sequence of N pulses, there are only N fixed 
points R^ in the path. Since between the pulses the Hamilto- 
nian H(t) = 0, we can define arbitrary orthonormal states as 
the basis \n' p ). Therefore, the geometric phase Ug so (T) is not 
unique. But notice that when we change the path (which is 
still fixed by the points R^), the error evolution U^uiT) also 
changes. A good choice of the path is the path that evolves 
smoothly along the points R^ and induces slow change of 
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FIG. 2: (Color online). The modulation functions F(if>) of a spin-^ in 
a scaled period. Red solid (blue dashed) lines for the real (imaginary) 
part of F((f>). (a) for the modulation function of a pulse sequence, (b) 
for a constant field B const , and (c) and (d) for fast varying fields B n 
and B 2n , respectively. 



Hf,%(R) and negligible error A£/e, t . From another point of 
view, we use R M to sample a given path, and when the number 
of sampling points is large enough the error AUeh is negligi- 
ble. The idea is illustrated by a two-level system in Fig.[TJ 

There will be extra errors in C/eitCO if there is an uncon- 
trollable Hamiltonian Ho + 0, besides the controllable fields. 
Hq may contain system-environment interaction, which is the 
case of an open quantum system framework. The system- 
environment interaction in general will induce decoherence 
in Ueh(T), for example, in conventional adiabatic evolution. 
However, if the dynamical decoupling control is implemented 
in a sufficient short time T, the effect of Hq can be negligible. 
Let U(T) be the evolution without Ho. The system evolution 

For 



-i J Q T U' l '(t)H U(t)dt 



with Ho is tWr) = U(T)T exp 

a bounded Ho, U tota j(T) — > U(T) when T — > 0. An interesting 
topic for future work is the design of dynamical decoupling 
control to suppress the effect of uncontrollable Hamiltonian 
more efficiently. 

A spin- j driven by a pulse sequence. An example of the 



pulses Eq. ( 16 1 is a sequence of ±n rotations along the direc- 

The 



tions x sin 6 cos 
parameter R^ = 



pf! + y sin & sin <^ + z cos & on a spin- 
(jpn). Each of the unitary pulse, 



2' 



P(<Pii) = ^ exp 



±j(-l)«W- 



(18) 



introduces a ±n phase shift between the instantaneous eigen- 
states |0+> = cos f e -*/ 2 | f ) + sin f e* /2 | |) and \(j>~) = 
-sin §e-* /2 | T) + cos fe* /2 | |). Here g(p.) = +1, and the 
dynamic phase is cancelled when the numbers of +n and -n 
pulses are equal. There are other methods to cancel the dy- 
namic phase using pulse sequences 126112711 . 

The geometric (Berry) phase from the evolution <f>o = 

to cf> T is U s ^(4>t) = 2± l0r><^ol e±,COSflfr/2 - And the error 
evolution 



''C 5 f W sine ( cicos# li><TI+h.c.)rf0 



(19) 
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FIG. 3: (Color online). The averaged error AUsn is suppressed by 
CP sequences with different pulse numbers N, with red circles (blue 
squares) for even (odd) N. Here (f> = 0, <f> T = 2n, and = n/6. 



where the modulation function F(<f>) = (— iy when (f> e 
(<f>p-i,4>n] [see Fig. |2ja)] . Note that if we apply 2n rotations 
on the spin-l, even though the energy gaps are larger during 
the control, the modulation function F(0) = 1 is a constant 
and the adiabatic evolution is not realized. 

In Fig.[3]we show the error AE/eh [Eq. ( 15 i] for the case of 
<pj — 2n and 9 — n/6 under the control of Carr-Purcell (CP) 
pulse sequences [28 ] 



<Pli = (<Pt - 0o) 



2fx- 1 
2N 



+ <po, for fi=l,. . .,N. 



(20) 



Here the CP pulses equally sample a circular path with the 
spacing 2n/N. Figure [3] shows that as the number of pulses 
increases, the error evolution is smaller, consistent with our 
theoretical result. The CP sequences with even number of 
pulses have better performance than the sequences with odd 
N. Note that when 9 — n/2, the error is completely eliminated 
(i.e., At/Err = 0) by the CP sequences regardless of the pulse 
numbers N. 

A spin-l^ driven by continuously varying fields. Fast 
varying fields that are changing continuously can also lead 
to adiabatic evolution and can have better performance 
than slowly varying fields. For example, the Hamilto- 
nian of a spin-j is H(t) = \o~ ■ B(t), where B(t) - 
B(t) [x sin 9 cos (f> + y sin 6 sin + z cos 9] with the parameter 
<p = ojt changing in time continuously. We consider three 
cases of the field amplitudes B(t) for comparison: (i) B n {t) = 
f [1 - ycos (QQ], (ii) B 2 „{t) = 2B n (t), and (iii) B const (f) = 
y(2 + y 2 )/80, which has the same average energy as B n {t) 

r— r — 

[i.e., J a |Sconstl 2 * = J a \B„\ 2 dt]. During the evolution there 
are degeneracy points when B„(t) = 0. The field B n {B2 n ) con- 
tributes a n (2n) phase shift in each period of 2n/Q.. We set 
y ~ 2.34213 so that the average of the modulation function 
e'X B M ds vanishes in a half period n/Cl. 

In Fig. |4] we plot AJ/eit for B„, B2 n , and B const as a func- 
tion of N' = Q./2n with the total evolution time T — 1 and 
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FIG. 4: (Color online). The errors At/gn is suppressed by continuous 
driving fields, with ipo = 0, i/>t = 2n, and = n/2. For 6^, Af' = 
Q/(2n) is the number of applied n pulses. The results are shown at 
integer numbers of N', with the red squares (black triangles) for the 
fast field B„ (Bin) and the green circles for the field B const of constant 
amplitude. For other values of N', the results are shown for S const by 
blue dots and for B H by red empty squares from N' = 1 to 6. 



a) = 2n. For B„, the integer values of N' is the number of 
accumulated n phases during the evolution. Increasing N' 
(e.g., increasing the field amplitudes) is equivalent to increas- 
ing the evolution time in adiabatic evolution. As shown in 
Fig.|4j the fast varying field B n realizes the adiabatic evolution 
Ueu(T) — > / even though the field amplitude changes rapidly 
and the energy gap can be zero during the evolution. The field 
B„ with even number of n phase shifts is much more efficient 
than the slowly varying field Z? CO nst used in conventional adi- 
abatic evolution, because the modulation function e'/ B * df is 
more efficient than e ! J B «« dt (see the modulation functions in 
Fig. [2]). Even though ^2^(0 has larger amplitude and energy 
than B n (t), it can not realize adiabatic evolutions since the av- 
erage of the modulation does not vanish, J e' k Bl -« ds dt + 0. 
Thus larger field amplitudes (or energy gaps) do not always 
lead to better adiabatic evolutions. An interesting topic for fu- 
ture work is to design more efficient control fields that have 
better performance than the example B„(t). 

Conclusions. We have developed a general formalism to 
quantify the error transition in adiabatic evolution. We have 
shown that the instantaneous eigenenergies and the degener- 
acy of the Hamiltonian can change rapidly during the adia- 
batic evolution. Instead of changing the parameters R contin- 
uously and adiabatically, we can use a pulse to sample a num- 
ber of points on the path to realize adiabatic evolution. When 
the parameters R change continuously, fast varying fields can 
lead to better adiabatic evolution. 
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